How to adapt broad-band gravitational-wave searches for r-modes 
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Up to now there has been no search for gravitational waves from the r-modes of neutron stars 
in spite of the theoretical interest in the subject. Several oddities of r-modes must be addressed to 
obtain an observational result: The gravitational radiation field is dominated by the mass current 
(gravitomagnetic) quadrupole rather than the usual mass quadrupole, and the consequent difference 
in polarization affects detection statistics and parameter estimation. To astrophysically interpret 
a detection or upper limit it is necessary to convert the wave amplitude to an r-mode amplitude. 
Also, it is helpful to know indirect limits on gravitational-wave emission to gauge the interest of 
various searches. Here I address these issues, thereby providing the ingredients to adapt broad-band 
searches for continuous gravitational waves to obtain r-mode results. I also show that searches of 
existing data can already have interesting sensitivities to r-modes. 
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I. INTRODUCTION 

In recent years the LIGO Scientific Collaboration and 
Virgo Collaboration have performed numerous searches 
for continuous gravitational-wave (GW) emission from 
rapidly rotating neutron stars. These include narrow- 
band searches for known pulsars [H-Q (i-e. on or near 
a timing solution), broad-band searches for the (non- 
pulsing) accreting neutron star in the low-mass x-ray 
binary Sco X-l [a 0, an d all-sky broad-band surveys 
for previously unknown neutron stars [3 |8HT2| . Also, a 
broad-band search for the non-pulsing central c omp act 
object in supernova remnant Cas A is underway [13J. 

Cas A is the youngest object targeted by a contin- 
uous GW search so far, and at 300 years old it may 
still be subject to the r-mode GW emission mecha- 
nism. The r-modes [l4| are rotation-dominated oscilla- 
tions which are driven unstable by GW emission [15l . Il6j 
(the "Chandrasekhar-Friedman-Schutz" version of the 
two-stream instability fl7l. [ill ) even in the presence of 
nuclear- matter viscosity 19|-[21 1 . This makes r-modes ob- 
servationally a very interesting possibility for GW emis- 
sion in newborn neutron stars [22j and accreting neu- 
tron stars in low- mass x-ray binaries [23|, [24|. Nonlin- 
ear hydrodynamic saturation of the r-mode instability 
is now thought (under most conditions) to occur at am- 
plitudes several orders of magnitude below those origi- 
nally guessed (25j . making unlikely the most optimistic 
scenario of extragalactically detectable highly chirping 
signals from young neutron stars (22j . However a low 
saturation amplitude would also keep r-modes active in 
neutron stars u p t o thousands of years old under the 
right conditions [25| . Therefore it is interesting to search 
for continuous nearly-periodic r-mode GW emission from 
Cas A and other very young pulsing or non-pulsing neu- 
tron stars. 

However up to now the continuous GW observational 



literature has assumed that GW emission is from "ellip- 
ticity," i.e. an £ = m = 2 mass quadrupole which rotates 
with the star. By contrast r-modes are dominated by 
I = m = 2 current quadrupole emission [19]. The mul- 
tipole structure of the radiation field affects the wave 
polarization and detector response to a signal. These in 
turn affect detection statistics and parameter estimation, 
most notably the procedure used to convert the detector 
response h(t) into an intrinsic GW strain amplitude ho 
(or an observational upper limit on it). Also for an r- 
mode the astrophysical interpretation of a measured ho 
or upper limit must lead to a mode amplitude rather than 
an ellipticity. (If the mode amplitude evolves in a com- 
plicated way, as predicted under some conditions (2(| , an 
observational upper limit would apply to a time-rms am- 
plitude.) Finally, indirect limits on GW emission used as 
milestones for search sensitivities and to evaluate poten- 
tial searches are different for r-modes than for ellipticity. 

In this paper I address these amplitude- and 
polarization-related issues for r-mode GW emission. 
This is sufficient for broad-band searches, i.e. those not 
targeting a known pulsar. I do not address two im- 
portant issues related to frequency and phase: The r- 
mode frequency as a function of neutron star spin fre- 
quency depends on the equation of state and relativistic 
effects (e.g. Frequency estimates are needed for 

GW searches targeted at known pulsars and will be ad- 
dressed elsewhere (28|. Also, in some scenarios [26] the r- 
modes may not maintain phase coherence over likely ob- 
serving times. This also will be addressed elsewhere [29| . 

The rest of this paper goes as follows: In Sec. HH I 
briefly review the basics of r-modes and the multipolc 
formalism of GW generation. In Sec. IIII1 I show how 
these issues relate to GW observations, translating the 
theory-oriented quantities in Sec. |TT]to those used in the 
observational literature. I present indirect limits on r- 
mode GW emission in Sec. IIV1 I summarize in Sec. [V] 
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Unless otherwise noted I use geometrized units, where 
Newton's gravitational constant and the speed of light 
are unity, and 1 M Q = 4.92 x 10~ 6 s. I also use "neutron 
star" in a broad sense which includes the possibility of 
quark matter, meson condensates, and other exotica. 



where a is a dimensionless amplitude, R is the stellar ra- 
dius, and the magnetic-parity vector spherical harmonic 
is defined in terms of scalar spherical harmonics Y g - Ja 
as 30] 
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II. MULTIPOLE STRUCTURE 

Here I describe ellipticity and r-mode GW emission in 
the context of the multipole formalism for GW genera- 
tion, particularly the canonical review by Thorne [3(| • 

The LIGO and Virgo continuous-wave observational 
literature characterizes the mass quadrupole moment of 
a neutron star in terms of an equatorial ellipticity 
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where N k is the unit vector pointing from the center of 
the star. The GW frequency f = 2tt/oj is identical to the 
mode frequency, which is related to the spin frequency 
by El 



£ + 1 



(using £ — m for the proper r-modes). 

The GW strain tensor can be expanded as [3C 
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where the z-axis is the rotation axis of the star and I a f, 
is the moment of inertia tensor. The latter is identical to 
the mass quadrupole tensor up to a trace, which is not 
important here. The x and y axes are chosen to co-rotate 
with the star so that e does not oscillate. It is convenient 
to relate e to the scalar mass multipoles which are 
given in a non-rotating frame and the Newtonian limit 
by 



16tt 



'(£ + !)(£ + 2) 



(2£+l)UV 2£{£-l) 



d 6 rpr l Y 



(2) 



1 °° / rl \ 1 



1=2 ra=- 



B2,£,m 



where the spin- 2 tensor spherical harmonics of electric 
(E) and magnetic (B) parity are defined [30] 
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with p representing the mass density, r the radial coor- 
dinate, and the standard spherical harmonic. For 
a real-valued density perturbation fixed with respect to 
a star rotating with angular velocity J7, the time depen- 
dence of these is I e ' m = I l ~ m * cx e ~ mnt . From e.g. 
Eq. (4.7a) of Ref. [13] we obtain the relation 



r2 : 2| 



(3) 



This can be checked by examining the GW luminosity: 
Inserting Eq. © into Eq. (4.16) of Ref. (30| obtains, e.g., 
Eq. (20) of Ref. d: 
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where uj is the angular frequency of the wave, equal to 
mil for a perturbation rotating with the star. (The 
above expressions assume a negligible contribution from 
all mass multipoles except I 2,2 — I 2, ~ 2 *, but there are ar- 
guments that I 2,1 also can be significant in some stars — 
see Jones [3l| for a new one and a summary of older 
ones.) 

In Newtonian gravity and the slow-rotation approx- 
imation, an r-mode is an Eulerian velocity perturba- 
tion [33 
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As well as the mass multipoles I e ' m we encounter the 
gravitomagnetic or current multipoles 
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The time derivatives in Eq. ([8]) suppress GW emission 
from higher-^ multipoles by powers of a characteristic 
velocity, and thus the £ = 2 multipoles tend to be dom- 
inant. For r-modes, the £ = 2 mode is also the most 
unstable and the least damped by viscosity (l9| . 

It is also useful to write the full GW luminosity includ- 
ing all multipoles [3(| 



~ 32tt ^ ^ 
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dt 



(12) 

where the angle brackets denote a time average. 

By using Eq. (jll[) and orthonormality relations from 
Thorne j30|, we see that an I — 2 r-mode produces a 
current quadrupole 22] 



g2.2 = 
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where J is a dimensionless functional of the neutron-star 
equation of state and M and R are the mass and radius 
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of the (unperturbed) star. Neglecting the r-mode den- 
sity perturbation, all multipoles other than S 2,2 vanish — 
including I 2 ' 2 and I 2 '~ 2 . 

The r-mode amplitude a which I use here is related 
to others as follows: The results of Arras et al. [25[ are 

expressed in terms of Ai — J/2, sometimes called 
c a in other papers, which is about 0.1a for fiducial neu- 
tron star parameters (see below). Bondarescu et al. (2(| 
express most of their results in terms of C Q , which re- 
moves the adiabatic change of amplitude as the star spins 
down. For the example mode triplet they use, the zero- 
viscosity parametric instability threshold is normalized 
to |C Q |o = 0.011 s^ 1 / 2 . This makes their C a amplitudes 
about 1.1 x 10 3 a^fi/8.4 x 10 3 s^ 1 or about 200a at a 
spin frequency of 100 Hz. 

In Eqs. |5| and ([7]) and the others derived from them, 
I have neglected corrections of relative order H, 2 /(ttp) 
where p is the mean density of the non-rotating star. 
These include the density perturbation associated with 
an r-mode, which contributes to the GW tensor and mass 
quadrupole also at relative order Vi 2 /(np). These correc- 
tions rise to the 10% level only for stars rotating faster 
than the fastest known pulsars [2l[ . In view of the uncer- 
tain factors of 2 or more in quantities such as the moment 
of inertia (see below) , these and the somewhat larger cor- 
rections due to relativistic gravity [33[ are not important. 
The exception is the correction to Eq. ([7]), which affects 
narrow-band searches for known pulsars and will be ad- 
dressed elsewhere [28| . 

III. DIRECT OBSERVATIONS 

Here I discuss the effects of r-mode multipole structure 
on aspects of direct GW observations, including polariza- 
tions, detection statistics, waveform parameter estima- 
tion (especially amplitude), and upper limit procedures. 

Equation © implies that GWs from r-modes will have 
different polarizations than GWs from ellipticity. To see 
this compare two GWs with equal luminosity (fT2)) . one 
dominated by a mass multipole I l ' m and one dominated 
by the corresponding current multipole S e ' m . The linear 
+ and x polarizations arc projected out of the GW strain 
tensor as 

h+ = h jk (p j p k -q>q k ) /2 = h jk e j £ ) 

= h ]k (p>q k + q>p k ) /2 = h jk e? k , (14) 

where the unit vectors p° and q J are orthogonal to each 
other and to N J . Using the identity qj — tj a bN a pb with 
Eq. (HSJ), we find that 

rpB2,t,m jk _ _rpE2,i,m jk 
1 jk e + — 1 jk e x j 

rpB2,e,m jk _ rpE2,l,ra jk /-, r\ 

1 jk e x — 1 jk e + ■ l 10 J 

For a fixed luminosity, the GW tensor hj k of the current- 
dominated GW is obtained from that of the mass- 



dominated signal by taking T- k ' ,m — > Tj k ' ' m . Substi- 
tuting Eq. flSJ) into Eq. (fT4"|) and combining with Eq. ([15]) , 
we see that this takes 

(h+,h x )->(-h x ,h+). (16) 

Since the + and x polarizations are positive and nega- 
tive parity eigenstates, and the difference between mass 
and current multipoles is parity, this was to be expected. 
Note that the GW power radiated at a given angle from 
the star's rotation axis is proportional to h 2 ^ + h\ — e.g. 
Eqs. (4.12) and (4.13) of Thorne [H— and thus is not 
affected. 

The strain response of a detector is the contraction 
of the GW tensor hj k with a tensor describing the de- 
tector's beam pattern. It can be written by combining 
Eqs. (10)-(11) and (20)-(22) of Ref. [34] as (making the 
time dependence explicit) 

h(t) = [a(t) cos 2tp + b(t)shx2ip]h + (t) 

+ [b(t) cos2ip - a(t) sin 2ip]h x (t), (17) 

where ip is the polarization angle, explained simply in 
footnote 4 of Ref. [l|, whose definition is equivalent to 
fixing p J and q J in Eq. (|14l) . (The precise forms of the 
modulation functions a and b arising from the rotation 
of the Earth-based detector are not needed here.) From 
this we see that Eq. (fTl))) is equivalent to taking 

ip -> ip + ir/A. (18) 

The transformation (|18|) lets us quickly examine the suit- 
ability for r-modes of data analysis methods developed 
for ellipticity. 

Many GW search methods assume a uniform Bayesian 
prior on ip and thus are not affected by the transforma- 
tion JTHJ). The ^-statistic of Jaranowski et al. [HJ is one, 
as can be seen by the lack of ip i n their Eq. (55) and its 
derivation. The ^-statistic was derived as a frequentist 
maximization over ip (and other angles); but when de- 
riving a Bayesian alternative K-statistic, Prix and Krish- 
nan [35j explicitly showed that T has an implicit uniform 
prior on ip. The K-statistic itself has an explicit uniform 
prior on ip. The multi- interferometer /--statistic [36[ 
performs a weighted sum of T over detectors and thus 
also has an implicit uniform prior in ip. The heterodyn- 
ing method used in most known-pulsar GW searches is 
Bayesian and usually uses an explicit uniform prior in ip, 
as in Eq. (15) of Dupuis and Woan [33] ■ Semi-coherent 
searches for unknown pulsars based on combining the raw 
power [1, [^, [ll[ or .F-statistics [13, [HI °f short stretches 
of data also effectively use implicit uniform priors on ip 
(see Abbott et al. @, especially the appendices, for de- 
tails of the former). The radiometer search for Sco X- 
1 was adapted from a stochastic background search 
and explicitly assumes no preferred polarization. 

Some recent known-pulsar searches 0, [1[ use values 
of ip from measurements of pulsar wind nebulae. The 
measurement and error estimates can be folded into a 
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Bayesian method [37| as non-uniform priors on ip, or the 
best ip value can be inserted into the ^-statistic to ob- 
tain the ^-statistic [HI]. In order to cover the possibility 
of r- modes, such known-pulsar searches will need to tar- 
get not only different GW frequencies from the ellipticity 
case but also ip — > ip + it /A. 

Parameter estimation of candidate signals is also af- 
fected by Eq. (fT5|) . The observational literature charac- 
terizes the continuous GW amplitude with an intrinsic 
strain amplitude defined in terms of ellipticity as 



h = r 1 uj 2 I zz e = r~ 



(19) 



where the second equality uses Eq. ([3]). This ho is the 
amplitude of the response of a hypothetical detector at 
either of the Earth's poles to a signal originating from 
a star over either pole whose rotation axis is parallel 
to that of the Earth. It is also simply related to the 
GW luminosity. The full parameter estimation problem 
is lengthy and I do not address it here, but I note the 
following simple and useful approximation. A detected 
signal will be integrated for much more than one day, 
since even semi-coherent search candidates will be fol- 
lowed up coherently. In this limit Jaranowski and Krolak 
[391 ] performed detailed simulations confirming that the 
beam-pattern modulation averages out and the signal is 
effectively replaced by an unmodulated sine wave with 
amplitude 



h ee = ho^/A + B cos(4V>), 



(20) 



where the lengthy expressions for the functions A and B 
are given in their Eqs. (36) et seq. Thus ip — > ip + 7r/4 
simply flips the sign of B, which for most most source 
parameters is a few percent correction to h e g. (It may 
seem strange that hg is affected; this is because we are 
no longer transforming ip at fixed luminosity but rather 
for fixed detector response.) 

An accurate measurement of ip could yield information 
on whether a signal is from ellipticity or an r-mode: For 
a known pulsar the ratio of GW frequency to spin fre- 
quency already distinguishes between mechanisms and 
thus ip is redundant. But if the GW signal comes from a 
pulsar wind nebula without a timed pulsar, and the GW- 
measured ip (assuming ellipticity) is inconsistent with the 
orientation of a jet or torus, this would indicate r-mode 
emission. Also, a GW signal coming from Sco X-l or an- 
other non-pulsing accreting neutron star could be com- 
pared with possible orientations from radio jets or x- 
ray reprocessing [4l[ (since the difference between these 
estimates is less than 7r/4). 

Upper limits on ho from the uniform-^ searches are 
obtained in terms of populations of software- injected sig- 
nals uniformly distributed in ip, and hence the upper limit 
procedures for ellipticity remain valid for r-modes. Up- 
per limits on ho from searches with a given ip determined 
e.g. by a pulsar wind nebula will need to use populations 
of injections taking Eq. (fT5)) into account to produce sep- 
arate limits for ellipticity and r-modes. 



To convert an estimated or upper-limit ho to an r-mode 
amplitude, it is convenient to write the GW luminosity 

as 

dE 1 



-u?r 2 hl 
dt 10 

and compare it to Eq. (fT2"[) , 
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~dt 
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if S 2 ' 2 is the only non- vanishing multipole. [Note that the 
numerical coefficient is 1/2 that of Eq. (j4} because the r- 
mode is traditionally defined as a complex perturbation 
while the ellipticity is defined as real, and e lult contains 
twice as much power as cos(wt) which lacks the imaginary 
part.] Equating ([21]) and (|22|) obtains 




-1, ,3 



(23) 



This equation allows us to convert ho to a for a fiducial 
value of MR 3 J, in the same way that h is converted 
to e for a fiducial value of I zz = 10 45 gem 2 . (In both 
cases what is really measured is a quadrupole — and, if the 
star's rotation axis is known, the parity which determines 
if it is a mass or current quadrupole.) Much of the theo- 
retical neutron star literature uses a polytropic equation 
of state with polytropic index n — 1, which (for Newto- 
nian gravity) fixes J rs 0.0164 [2^ . Combined with the 
above choice of I zz and the usual choice of M = 1.4 Mq, 
this fixes R ~ 11.7 km rather than the common choices of 
10 km and 12.53 km. Inverting Eq. (|2"31 and substituting 
the fiducial neutron star values above, we obtain 



a = 0.028 



10- 



1 kpc 



100 Hz 
~1~ 



(24) 



Note that modern investigations of mode saturation [26( 
indicate that a may vary substantially over a typical ob- 
serving time. Since detection statistics respond to inte- 
grated signal power, the a inferred from them is really a 
time-rms average. 

Here and below, when I give numerical values I assume 
the fiducial neutron-star parameters above and I do not 
write scalings with M, R, I zz , and J. That is because, 
even if the equation of state (EOS) is not specified, as 
long as it is barotropic only two of the quantities are in- 
dependent. By noting that I zz can be written as MR 2 
times a relatively EOS-independent function (e.g. Lat- 
timer and Prakash [421 ]) except for very low- mass stars, 
and that J is less EOS-dependent than MR 3 , I estimate 
that the EOS-related uncertainty in these quantities is 
dominated by uncertainties in M and R and is therefore 
usually a factor 2-3. 



IV. INDIRECT LIMITS 

Indirect limits on GW emission inferred from electro- 
magnetic astronomical observables are useful to gauge 
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the astrophysical interest of existing and planned GW 
searches. As in Ref. [43{, where I summarize this and 
related issues in more detail, I divide the sources into the 
four categories of accreting neutron stars, known pulsars 
and non-pulsing stars without accretion, and the large 
population of unseen neutron stars sought by all-sky sur- 
veys. 

For known pulsars the primary indirect limit is the 
"spin-down limit," obtained by assuming that all of the 
observed change in spin frequency is due to GW emission. 
Thus Eq. (|2ip is equated to the kinetic energy loss 



— = i zz iiii 

dt 



(9/16)J„ww 



to obtain 



1 45/ zz P 



(25) 



(26) 



where P = 2tt/H, is the observed pulse period. For r- 
modes this is 3/2 times the value for ellipticity-dominated 
emission, given for instance by Eq. (2) of Wette et al. [l3[ . 
Substituting Eq. (|23|) obtains 
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1.4 



P 



5/2 



P 



1/2 



= 0.033 



10 ms/ \io- 10 

/100 Hz x 7/2 



V / 




where the last equality uses / = 4/3/ P. Again these 
are time-rms a limits if a fluctuates over the observation 
time as possible in many young neutron star scenarios 
[2(1 . Based on numbers from the on-line version 1.40 
of the ATNF catalog [44| , the pulsars with dipole spin- 
down ages P/(2P) less than 10 4 yr generally have a s d 
of order unity, except for J0537— 6910 which has about 
0.1. Reisenegger and Bonacic [45[ also pointed out that, 
under the right conditions of cooling and viscosity, r- 
modes could remain active in millisecond pulsars for a 
long time (comparable to the spin-down age) after the 
star stops accreting as a low-mass x-ray binary. The mil- 
lisecond pulsars typically have a s d of order 10 -6 -10 -4 , 
with a few as high as order 10 -2 . 

A numerically stricter but more model-dependent indi- 
rect limit was obtained by Palomba [461 ] for a few known 
pulsars by incorporating additional information on the 
age and braking index 



(28) 



If £1 is proportional to a power of f2, the braking in- 
dex picks out that power. For radiative braking from a 
static (electromagnetic or GW) multipole £, the power 
is n = 11 + 1, while for particle winds it can be lower. 



The highest value usually considered is n = 7, which cor- 
responds to constant a. As first argued from adiabatic 
invariance |47j and later seen explicitly in nonlinear hy- 
drodynamic saturation calculations [26j . the long-term 
average of a might scale as SI" 1 / 2 , resulting in n = 6. 
Palomba [46| took to be a sum of two powers of fl, 
one power 5 (GWs from ellipticity) and one a free pa- 
rameter; and performed numerical spin-down evolutions 
of the pulsars over a wide parameter space constrained to 
be consistent with their known ages and present values 
of f2, Cl, and n. Because all observed braking indices are 
less than 3, the f2 5 GW component of fl is constrained 
to be less than the spin-down limit by some factor — in 
the case of the Crab pulsar, this limit on ho is 2.5 times 
stricter than the spin-down limit. Without performing 
such detailed simulations, it is clear that since r-modes 
have fl proportional to higher powers of O, this type of 
limit on the Crab would be stricter than the spin-down 
limit by more than the factor 2.5 for ellipticity. 

Compare a s( j to the theory of r-mode non-linear hy- 
drodynamical saturation: The lowest zero- viscosity para- 
metric instability threshold (C a = 1 in the notation of 
Bondarescu et al. |26jj) corresponds to a a few times 10 -3 
for the frequencies of interest. The real threshold, which 
tends to serve as an attractor for r-mode evolutions, de- 
pends on temperature as well as frequency and can be 
an order of magnitude higher (e.g. their Fig. 6). But it 
is still below the values of a s d for young pulsars, which 
are high enough to lie in the "run-away" regime of non- 
linear hydrodynamics which requires further study [26| . 
For millisecond pulsars, the appropriate comparison is to 
lower values (see the discussion of accreting stars below) . 

For non-accreting neutron stars without observed 
spins, such as the central compact object in Cas A, age- 
based indirect limits are obtained by substituting the age 
of the object into the spin-down limits under the assump- 
tion that the object has spun down significantly and pre- 
dominantly by GW emission [131 ] . In this case we have 
the relation 



P/P ={n- l)o, 



(29) 



where a is the age of the neutron star and the braking 
index n is assumed to be constant or appropriately aver- 
aged. In parallel to the derivation of the ellipticity ver- 
sion of this limit [11], we substitute Eq. into Eq. (j2"6]l 
to obtain y/Z/2 the value for ellipticity, or 



C = 1.5xl0- 24 



1/2 



3.4 kpc 



6 



1 



1/2 



(30) 



(The fiducial values are for Cas A and constant-a evo- 
lution.) Similarly for the age-based indirect limit on r- 
mode amplitude we obtain 



15 J^I^M-'R- 3 



80? 



(n-1) 



-1/2 



= 0.14 



300 yr 



1/2 



100 Hz 
~1~ 



1/2 



(3T) 
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where again the limit is a time-rms value if a fluctu- 
ates quickly compared to the spin-down timescale. For 
Cas A, the youngest known neutron star, the values of 
a age are 0.14-0.005 over the 100-300 Hz band considered 
by Wette et al. 13]. Several similar objects, e.g. the non- 
pulsing ones listed by Halpern and Gotthelf [Hj, have 
indirect limits of the same order of magnitude. At a few 
hundred Hz these values are comparable to the paramet- 
ric instability thresholds in several scenarios analyzed by 
Bondarescu et al. [26[ , with zero- viscosity a thresholds a 
few times 10 -3 and true thresholds several times higher. 

For accreting neutron stars the standard indirect limit 
is derived from the argument that accretion torque and 
GW torque are in equilibrium, originally made by Pa- 
paloizou and Pringle [49j , put on firm astrophysical foot- 
ing by Wagoner [50j, and later tied to observations of 
spins of accreting neutron stars [23|, [24| . From dJ/ dt = 
(2/ui)dE/dt—as in Eq. (4.23) of Thorne we have 



dJ/dt, 



(32) 



equivalent to Eq. (4) of Watts et al. [51]. The indirect 
limit on ho is obtained by assuming accretion of Keple- 
rian angular momentum at the stellar surface and 100% 
conversion of gravitational potential energy to x-ray flux 
Ft measured at Earth to obtain 



i^acc 
"0 



/20^ CJ -l/2 F V2 M -l/4 jR 3/4 



3.3 x 10" 26 



3.9 x 10~ 7 erg cm -2 
/800 Hz x 1/2 



1/2 



(33) 



Here the numerical value is scaled to the average bolo- 
metric flux of Sco X-l (which does not pulse) and the 
highest spin rate observed in an accreting millisecond 
pulsar, both from Watts et al. [51]. For a given GW 
frequency /i§ cc is the same as the indirect limit for el- 
lipticity GW emission, but it is different for a fixed spin 
period. Combining with Eq. (|23[) we obtain 

135^ 



2048tt 7 /2 
5.1 x 10~ 6 



Fl^P^rM-^R-^j- 1 



F,. 



2.8 kpc 



3.9 x lO- 7 ergcm- 2 
800 Hz \ 7/2 



1/2 



(34) 



For comparison, Bondarescu et al. 52] find that r-mode 
evolutions of accreting neutron stars have a ps 10~ 4 
if there is no run-away (here the viscosities and thus 
the parametric instability thresholds differ from those of 
young neutron stars). That nominally corresponds to 
P = 4 ms or / ps 300 Hz for Sco X-l at the indirect 
limit. Since the r-mode amplitudes in these saturation 
studies are uncertain by at least a factor of a few due to 
uncertainties in stellar structure and damping rates, this 



is not a precise prediction of the spin period but rather a 
statement that Sco X-l is consistent with having a short 
spin period regulated by r-modes at or near equilibrium. 
(This was independently found in a more detailed ex- 
amination of possible r-mode evolutions J53].) However 
non-linear mode evolutions may avoid accretion torque 
equilibrium altogether, may spin up as well as down, and 
may go through intermittent episodes of r-mode activity; 
and realistic accretion may be more complicated than 
what is usually assumed. Therefore indirect limits from 
accretion equilibrium are much softer than the spin-down 
and age-based limits, which are derived from energy con- 
servation. (This is a good thing, since the former are 
much more pessimistic than the latter (43J.) 

All-sky surveys for continuous GWs are subject to sta- 
tistically estimated indirect limits derived from assump- 
tions about the galactic supernova rate and distribution 
of neutron-star parameters at birth. A simple version 
due originally to Blandford (unpublished) is described in 
Abbott et al. @ and is not changed much by r-modes 
since it essentially relies on /iq S ° and a planar galactic 
geometry to yield a population estimate of the brightest 
signal /i0° p a few times 10 -24 independent of frequency 
and the poorly-known ellipticity. Simulations by Knis- 
pel and Allen [54| using realistic spatial distributions of 
neutron stars produce ti^° v values which do depend on 
ellipticity and frequency. Investigating the effect of r- 
mode emission on this type of indirect limit estimate is 
a substantial undertaking which I do not attempt here. 
However r-mode emission has a steeper frequency de- 
pendence than ellipticity emission, resulting in r-modes 
dominating dE/dt for 



a > 



3V5 I zz 
16tt MQR 3 f 



87e 



P 



10 ms 



(35) 



(obtained by equating I 2 ' 2 to S 2 ' 2 ). This condition is sat- 
isfied even by the lowest parametric instability threshold 
for the frequencies and ellipticities contributing to the 
main result (Fig. 5) of Knispel and Allen 54[. Thus the 



presence in a similar simulation of a significant popula- 
tion with active r-modes should spin down stars more 
quickly into the best frequencies for GW detection and 
to produce a higher /iq OP , even assuming the scenarios of 
[26( and not invoking a low- viscosity run-away. 

The last form of the indirect limit (|27p for known pul- 
sars can also be used to interpret direct upper limits from 
all-sky surveys by overlaying the contours of a s d(/, /) on 
a plot similar to Fig. 41 of Abbott et al. [9j. Following the 
discussion around that figure we can say: The range of 
the multi-interferometer Hough search from Abbott et al. 
[|] was about 1 kpc at / = 100 Hz and |/| = 10~ 8 Hz/s, if 
/ for a star is GW-dominated — and Eq. (|2"7|) tells us that 
would require a ps 0.3. A star at / = 1 kHz with the same 
/ would have been detectable up to 100 pc (comparable 
to the closest known neutron stars) if GW-dominated, 
and this would require a ~ 1 x 10~ 4 , which is well within 
the range of possibilities for non-linear evolutions. 
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V. DISCUSSION 

I have shown that current data analysis methods can 
detect or set valid upper limits on continuous GW emis- 
sion from r-modes, in some cases with small modifica- 
tions. I have derived relations needed to interpret GW 
observational results in terms of r-mode emission. This 
allows broad-band searches for continuous GWs to infer 
r-mode amplitudes or upper limits. 

Searches for r-modes in known pulsars will also re- 
quire information on the ratio of GW frequency to spin 
frequency [HI], and will be more sensitive to the issue of 
r-mode phase coherence time due to their long integra- 
tion times [29(. 

I have also re-derived some commonly-used indirect 
limits on GW emission for the case of r-modes. Several 
young pulsars and non-pulsing neutron stars are inter- 
esting targets for searches for r-mode GW emission, and 
all-sky surveys can have interesting ranges, even with 



presently available LIGO and Virgo data. More young 
objects further away will be interesting with advanced 
LIGO and Virgo. 
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